Electronic states, Mott localization, electron-lattice coupling, and dimerization for 

correlated one-dimensional systems. II. 
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We discuss physical properties of strongly correlated electron states for a linear chain obtained 
with the help of the recently proposed new method combining the exact diagonalization in the Fock 
space with an ab initio readjustment of the single-particle orbitals in the correlated state. The 
method extends the current discussion of the correlated states since the properties are obtained 
with varying lattice spacing. The finite system of A'^ atoms evolves with the increasing interatomic 
distance from a Fermi-liquid-like state into the Mott insulator. The criteria of the localization are 
discussed in detail since the results are already convergent for N > 8. During this process the Fermi- 
Dirac distribution gets smeared out, the effective band mass increases by 50%, and the spin-spin 
correlation functions reduce to those for the Heisenberg antiferromagnet. Values of the microscopic 
parameters such as the hopping and the kinetic-exchange integrals, as well as the magnitude of both 
intra- and inter-atomic Coulomb and exchange interactions are calculated. We also determine the 
values of various local electron-lattice couplings and show that they are comparable to the kinetic 
exchange contribution in the strong-correlation limit. The magnitudes of the dimerization and the 
zero-point motion are also discussed. Our results provide a canonical example of a tractable strongly 
correlated system with a precise, first-principle description as a function of interatomic distance of 
a model system involving all hopping integrals, all pair-site interactions, and the exact one-band 
Wannier functions. 
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PACS Nos. 71.10.Hf, 71.27.+a, 71.30.th 
I. INTRODUCTION 



In spite of the enormous successes of the apprjOach 
based on the effective single-particle wave equationEJ for 
many 3-dimensional metals and semiconductors, the un- 
derstanding of the so-called correlated fermionic systems 
is still lacking, particularly for the systems of lower di- 
mensionality, d = 1 and 2. This is because in their de- 
scription of the electronic states the role of the long-range 
Coulomb interaction [is crucial, as the charge screening 
becomes less effectivaa. In result, the interaction cannot 
be regarded as small and phenomena such as the spiru 
charge separation^, the Mott metal-insulator transitioncl 
or the strong electron-lattice iuteraction (leading among 
others to the Peierls distortedu state) appear in various 
quasi-one-dimensional systemS|-| In addition, a normal 
metallic state in quantum wiresQ and nanotubeaD, as well 
as the superconducting state in organic Id metals are also 
observed, raising the question about the relative role of 
the single-particle dynamics and the repulsive Coulomb 
interactions treated in a nonperturbational manner. To 
phrase it differently, in one-dimensional systems we con- 
sider here the single-particle states and the Coulomb in- 
teractions are treated on equaLfooting from the outset. 

We have recently proposedBcl (those References are 
regarded as Part I) a new approach to the correlated 
fermion systems, which provides rigorous results for 
model (one-band) systems containing iV < 12 atoms. 
The results concerning the principal local characteris- 
tics are very rapidly converging and provide a reliable 
estimate of the system characteristics such as the micro- 



scopic interaction parameters, the ground state energy, 
the magnitude of Peierls distorsion or the zero point 
motion, all as function of the interatomic distance. It 
is the purpose of this paper to provide a complete pic- 
ture of the electronic and lattice properties by providing 
the band effective masses, the magnitude of exchange in- 
teractions and related to them the spin-spin correlation 
functions, the detailed discussion of a gradual transition 
from a metal to the Mott insulator, and the constants 
of the local coupling of electrons to the lattice. We also 
consider the dimerized state of the system. We believe 
that analysis of this type containing a combination of a 
rigorous treatment of the interactions in the Fock space 
combined with ab initio treatment of the single-particle 
wave-functions in the Hilbert state, that are allowed to 
relax in the correlated state, provides the first step in a 
rigorous quantitive analysis of linear-chain of finite size 
and correlated quantum dot systems. In other words, 
our first-principle method is specifically designed to treat 
such correlated systems, although some general conclu- 
sions are also discussed. 

The structure of the paper is as follows. In the next 
section (and in Appendix A) we summarize our method 
and provide the details not published before. In Sec- 
tion III the crossover from metallic to the Mott insulat- 



ing state is dealt with to provide the proper ground state 
for a subsequent analysis. In Section |^ we discuss the 
magnetic (kinetic-exchange) interactions and the spin- 
spin correlation functions. In Section ^ we determine the 
local electron-phonon coupling constants and compare 
their magnitude with that of the magnetic interactions. 
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The main feature of the above analysis is to provide the 
properties as a function of the lattice spacing. In this ■ 
spect our approach differs from the numerous solutionsc 
of parametrized models such as the (extended) Hubbard 
model, where the physical properties are discussed as a 
function of the model parameters. Here all one- and two- 
site parameters are calculated explicitly and this feature 
allows for a consideration of e.g. one-band system with 
inclusion of both the exact Wannier functions, all hop- 
ping integrals and all two-site interactions. 



II. METHOD 



The quantity a, which represents the inverse orbital size 
(and calculated in units of the Bohr radius ag) will 
be determined by optimizing the orthonormal atomic 
(Wannier) basis in the correlated state. In effect, the 
parameters are ea ~ {wi\T \wi), ti-j — {wi\T \wj), 

U = {WiW,\ Vi2 Ki^j = {w^Wjl Vi2 \w,Wj), Ji-.j 

{wiWj\Vi2\wjWi) = {wiWi\Vi2\wjWj), and Vi--j 
{wiWi\Vi2\wiWj) . The operator T represents the full 
single-particle lattice potential, i.e. 



r(r) 



2m 



N-l 

E 

]=0 



We start with the Hamiltonian containing N lattice 
sites with all hopping integrals tij and with all two- 
site interactions, which for the linear chain with periodic 
boundary conditions can be written as 



JV-l 



i-1 

E 

3=0 



E 



{Uj +V^j{nis+njs)) {al 



h.c. 



1} 



(1) 



The first term represents the atomic energy (we in- 
clude it explicitly, since changes with the varying lat- 
tice constant). The second describes the intraatomic 
Coulomb interaction (the Hubbard term). The next two 
terms represent the direct intersite Coulomb interaction 
Kij) and the Heisenberg-Dirac exchange (~ Jij). The 
fifth and the sixth term express respectively the single- 
particle ('~ tij) and the correlated- hopping (~ Vij) terms, 
whereas the last includes the pair-electron hopping i ^ j. 
The microscopic parameters have been defined before (cf. 
Appendix B of Ref.B) . Their values as a function of the 
interatomic distance R have been determined beforeO; 
here we provide in Table ^ the values of the interaction 
parameters together with their asymptotic analytic ex- 
pressions, which reproduce well their values for R^Agq 
(>2A). 

One should underline that the parameters are here de- 
fined in terms of exact Wannier function for this one-band 
(one-orbital-per-site) system, that is defined as 



N 



.(r) = ^A_,vE',(r), 



(2) 



where 5'j(r) represents atomic function of s-type cen- 
tered on site j: 



*,(r)^*(r~R,) 



3\ 1/2 



exp (— a |r — Rj l) . (3) 



E- 



R,7 



(4) 



where a.u. means the expression in atomic units. V12 — 
e^/|ri — is the usual Coulomb potential (we do not 
include any screening by e.g. core electrons as we want to 
discuss the model situation, but in a rigorous manner). 

TABLE I. Microscopic parameters values for various inter- 
action terms in the starting Hamiltonian for N = 8 sites. 



R [ao] 


U [Ry] 


Ki [Ry] 


K2 [Ry] 


Kz [Ry] 


2.0 


2.301 


1.077 


0.676 


0.450 


2.5 


1.949 


0.843 


0.499 


0.331 


3.0 


1.717 


0.692 


0.391 


0.259 


4.0 


1.452 


0.508 


0.269 


0.179 


5.0 


1.327 


0.403 


0.206 


0.138 


R 00 


1.250 


2/R 


l/R 


2/3i? 



R [ao] 


Ji [mRy] 


Vi [mRy] 


V2 [mRy] 


2.0 


9.54 


-18.07 


33.58 


2.5 


7.39 


-17.45 


19.58 


3.0 


5.59 


-16.08 


11.95 


4.0 


2.90 


-12.92 


4.49 


5.0 


1.26 


-9.64 


1.56 


R ^ 00 


'-^e-^logR 


2Re-^ 


4i?e-2« 



The inclusion of the full Coulomb potential (^) will in- 
troduce three-site terms already in the hopping integral 
t[j in the atomic basis; this feature follows the definitions 
and expressions for all the interaction parameters for the 
Wannier function (^) provided in Appendix A (RefiB does 
not contain all parameters explicitly and the evaluation of 
tij is absent, see Appendix B). Also, the outlines of the 
diagonalization method and of the single-particle-basis 
optimization are briefly summarized in Appendix C. In 
what follows we concentrate on the physical properties 
of the results for a linear chain with periodic boundary 
conditions containing up to iV = 12 atoms. Particular 
emphasis is put on the asymptotic properties (i.e. those 
weakly dependent on N) of those truly nanoscopic sys- 
tems depicted schematically in Fig. |l|. 
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III. METALLIC AND MOTT LOCALIZED 
STATES 

A. Band energy and the effective band mass 



We begin with the whole analysis of the band elec- 
tron properties. The values of atomic level position 

— J d^rw*{r)T{r)wi{r) and of the hopping matrix el- 
ements tm = J (i^ri£*(r)r(r)u'i+m(r) with m = 1, ... ,5 
are listed in Table || for R/ao = 2^5. At most, three 
first hopping matrix elements are important; this is suf- 
ficient to achieve the asymptotic properties already for 
N > 8. Such a fast diminution of tm with increasing 
m is due to the reduction of the spread of the Wannier 
functions {|Wi(r)} due to the intersite repulsive Coulomb 
interactionO. The space profile of the Wannier function 
along the chain direction for different interatomic dis- 
tances is shown in Fig. ^. The wave-function renormal- 
ization is apparent for smaller R/aQ, where metallic state 
is expected to appear. 

From Table || one can also see that the atomic en- 
ergy is strongly dependent on the distance R and even 
for R = 5ao it is decisively lower (ss —2.7 Ry) than that 
in the Is state of hydrogen atom (—1 Ry). This decrease 
is caused both by the long-range nature of the attractive 
potential V{r). 

The band energy can be calculated directly by intro- 
ducing the corresponding single-particle energy 



Af-l 
m=0 



/ 2Tnnkr, 



(5) 



where fc„ = -[N/2], ...,[{N- l)/2] {[x] denotes the in- 
teger part of x) is the quantum number (related to the 
wavevector via relation k — 2TTkn/NR) in the first Bril- 
louin zone. The profile of the bands evolving with the in- 
creasing lattice parameter R is shown in Fig. ^b. Only 
for the distances R^iaQ we can approximate the band 
with one hopping parameter ti, when the bandwidth is 
strongly reduced, as one would expect on the basis of 
the tight- binding-approximation (TEA) . The Fermi level 
is always at the point kR/n = 0.5, as we have assumed 
that we have one electron per atom. 

From the band energies we can obtain the effective 
band mass at either the band center {k — 0) or at the 
Fermi point {kp = ■n/2R): 



dk^ 



— h^kf 



k=0 



k—kf 



2 

i?2 



fe=0 



(6) 




k=ixl2j 



The relative efi'ective mass m* /me (where me is the 
free rest mass) are plotted in Fig. ^b (the inaccura- 
cies are due to the numerical differentiation of 6^). The 
mass grows with the increasing interatomic distance and 
reaches about 50—70% higher value (then me) for Rja^ « 
5. One should also note that even thought the Wannier 
functions are obtained from optimizing the energy of the 
interacting state, the masses are light (m* < me) for 
R < 3.5ao. What is more important, they practically do 
not depend on N for TV > 6. Therefore, in Fig. ||we have 
plotted the profile of ni* /me only for = 10 versus both 
R/ao and kR/n. 

The calculated band energies will serve as an input in 
the discussion of the relative role of the Coulomb inter- 
action at the onset of the Mott localized state. This is 
considered next by determining first the bare bandwidth 

W = 2 X]m=i and comparing it with the magnitude 
of the effective Coulomb interaction. 



B. Onset of the Mott localized state 

The determination of the basis {?i'i(r)}i=i....^jv allowed 
for determination of the interaction parameters (cf. Ta- 
ble |) and the hopping integrals tm (cf. Table ||) . There- 
fore, the bare bandwidth can be also defined as 

ly = max{efc} - minjefe}. (8) 

In Table III we present (for = 8) the effective orbital 
size ttH = ctmln i^^ units of ao), the bandwidth W, the 
effective interaction parameters U — Ki, the W/ {U — Ki) 
ratio, and the product of the carrier density nc — 1/R 
and the optimal orbital radius an = Q^min- pr<M" 
uct ncaH illustrates the Mott criterion for localizationll3, 
which can be generalized first to the case of d-dimensional 
infinite lattice, which takes the form n^i'^aH ~ 0.2. The 
critical value of the product 0.2 is reached for R k, 4.5ao 
and hence we should ask, whether this is a coincidence or 
if it reflects the localization onset for the correlated elec- 
trons in these nanoscopic systems. This criterion is not 
strongly dependent on N , as shown in Fig. || and hence 
represents an intrinsic property, independent of the sys- 
tem size for > 8. 



TABLE II. Single-particle parameters versus interatomic 
distance for A = 10 sites. 



R/ao 


fa 


ti 


t2 


is 


t4 


is 




[Ry] 


[Ry] 


[mRy] 


[mRy] 


[mRy] 


[mRy] 


2.0 


-4.461 


-0.585 


87.0 


-8.93 


1.29 


-0.413 


2.5 


-4.077 


-0.330 


44.0 


-4.10 


0.54 


-0.154 


3.0 


-3.712 


-0.200 


23.6 


-2.00 


0.23 


-0.006 


3.5 


-3.399 


-0.127 


13.0 


-0.99 


0.10 


-0.002 


4.0 


-3.138 


-0.083 


7.5 


-0.54 


0.05 


-0.009 


4.5 


-2.920 


-0.055 


8.2 


-0.66 


0.05 


-0.008 


5.0 


-2.737 


-0.037 


4.3 


-0.28 


0.02 


-0.002 



(7) 
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For that purpose we have calculated basic quantities 
signalling such a crossover, each of which can be charac- 
terized briefly under the following headings: 

(i) The total spin length per site. We characterize it 
by its square, i.e. by (S^) = (0| |0), where = 

(5+,S'r,S'f) = {a\^a^l,a\^a^^,,{n,^^ -n,^)l2^ is 

the electron spip-,operator for the atomic site i. It 
is easy to provetll that 



(9) 



Therefore, in the atomic limit with one electron 
per atom we have that <P = (ni|ni|) — and 
hence the spin length acquires the Pauli-spin limit 
S{S + 1) = 3/4. In the opposite, Hartree-Fock 
limit (C/eff < W), (n,Tn,a) = 1/4, and (Sf) = 3/8. 
Hence, the quantity 0m = (4/3) (S|) starts from 
the value 1/2 for i? and approaches the value of 
unity when the atomic-like localized states are more 
proper. Note that the ground state of the system 
is always a total spin singlet, i.e. (0| Ylif=i I) — 0- 

(ii) The dispersion of the statistical distribution of 
function nua = (0| a\^ai„ |0} = (0| hka |0). The dis- 
tributions in both the optimized correlated state 
(i.e. with a = amin) and in the state without such 



an optimization (i.e. for a 



) is displayed 



in Fig. 0ab. As the case with the optimized size 
= a~?^^ has a lower energy (see below), the be- 
havior of confirms the existence of the Fermi 
ridge aX k — kp — n/2R for R ~ 2ao followed by 
its gradual diminution to its disappearance with 
the increasing R. The presence of the Fermi ridge 
speaks directly in favour of delocalized (metallic) 
state of electrons0. The distribution for i?^4ao 
is essentially the Fermi-Dirac distribution modified 
by the Fermi liquid effects, which is smeared out 
completely for R > Soq. To put this process on 
the quantitive grounds we have calculated first the 
dispersion of the statistical distribution 

k<7 



2N ^ 



I flka |0) 



(10) 



In the Hartree-Fock limit this quantity can be cal- 
culated by assuming that {fika) = 0(/^ — ffe), so 
that u'^inka} = 1/4, whereas it vanishes for the 
even momentum distribution {uka = 1/2), when 
the particle position is sharply defined on atom, 
i.e. for the localized states reducing to the atomic 
states. 



(iii) The nearest neighbor spin-spin correlation func- 
tion &AF — — i^i ■ Si+i). It should be zero in the 
Hartree-Fock limit and reach the value (3/4) for the 
nearest-neighbor singlet state which mimics the an- 
tiparallel orientation of the classical spins. 

(iv) The Fermi discontinuity (ridge) disappearance. It 
is defined ag|l3 



nk=kf 



(11) 



This difference has been interpolated by a parabola 
Uka = Oik"^ -\- Pk -\- 7 from both sides leading to 
the value An^^"^ at the jump. Such an interpo- 
lation simulates a quasicontinuous function Uka, 
which would be present for large N. The resul- 
tant R dependences of both computed and inter- 
polated Ankp values for iV = 10 are listed in Ta- 
ble IV. The discontinuity disappears in the range 
of i?c = 4 4.5ao. The large uncertainty is due to 
the poor statistics of the points (3 points on each 
side). However, it is close to the values deducted 
from R dependence of the quantities (|l)-([i|), as we 
discussed next. 

The quantities characterized (p|)-(pi|) are displayed in 
Fig |[ where we have shadowed the difference between 
the results for = 6 and 10 (the results for N = 8 fall 
in the area) to amplify the convergence of the numeri- 
cal results. We see that the quantities 6mi, ©af, and 
Qm acquire the atomic values within the 5% range for 
R/ttQ « 5, which corresponds to the interatomic distance 
R « 2.6A. It should be underlined that we do not ex- 
pect any discontinuous phase transition for this finite (in 
fact, nanoscopic!) system. However, with the help of 
the characteristics provided above one can define an ex- 
perimental criterion of localization (here we suggest that 
the 5% margin for the characteristics to fall within the 
atomic limit values, is a natural one). Also, all the char- 
acteristics defined above are defined from the side of the 
delocalized (metallic) state. To define properly the Mott 
insulating state as that of a Heisenberg magnet (i.e. with 
frozen orbital degrees of freedom), we have to consider 
the spin-spin correlation function directly, as well as the 
magnitude of the superexchange (kinetic exchange) , both 
as a function of interatomic distance. 

TABLE III. Effective Bohr radius {an), the bare bandwith 
(W), the interaction parameter {U — Ki), the bandwith to in- 
teraction ratio, and the Mott criterion (ncan), respectively. 



R/ao 


an/ao 


W 


U-Ki 


W/{U-Kj) 


ncan 


2.0 


0.570 


2.381 


1.224 


1.945 


0.285 


2.5 


0.667 


1.340 


1.106 


1.212 


0.267 


3.0 


0.750 


0.810 


1.026 


0.790 


0.250 


3.5 


0.818 


0.512 


0.976 


0.525 


0.234 


4.0 


0.874 


0.333 


0.944 


0.353 


0.219 


5.0 


0.948 


0.148 


0.925 


0.160 


0.190 
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The Fermi ridge is sharply defined in the many-fermion 
system only when thje perturbation expansion is conver- 
gent (cf. Luttingerll3) . It is certainly not at the Mott 
metal-insulator transition. Our results show that the 
Fermi discontinuity disappears with the increasing lattice 
parameter even in the situation when we have a crossover 
transition from a metal to an insulator. This means that 
low-dimensional (finite) systems cannot always be ana- 
lyzed perturbationally even thought they do not exhibit 
phase transformation in the thermodynamic sense. Ob- 
viously, part of the jump at the Fermi points is due to 
the discretness of the particle quasimomentum in this 
finite-size system. 



IV. KINETIC EXCHANGE AND SPIN-SPIN 
CORRELATIONS 

In our system the total spin is conserved. Since 
the ground state is a spin singlet, we have that 



|0) = 0. From this we can derive the sum 

rule for the ground state of the Heisenberg system in the 
form 



JV-l 

E 

m— 1 



(Si • Si^^ri) — —Slot (5'tot 



1) = 0. (12) 



In our case the spins are defined in the Fock space (see 
above). Therefore, the sum rule may not be obeyed, 
since the spin magnitude is smaller. This is explicitly evi- 
dent in Table where the spin-spin correlation functions 
{Si ■ Si+m) are listed as a function of R. The long-range 
correlations set in and oscillate in sign as the atomic limit 
is approached . This feature of the correlation function 
must be induced by the Anderson antiferromagnetic ki- 
netic exchange interaction. We have calculated the ki- 
netic exchange integrals between m-th neighbors defined 
aM 



J, 



(m 
kex 



) _ 4 (t^ + Vm) 



(13) 



and displayed them in Table VI. The total strength 
of kinetic exchange {~ Yl^m=i ^vl'^) ^^^'^ been com- 
pared with the corresponding quantity (« X]m=i "^m) 
pressing the strength of the Heisenberg-Dirac exchange 
(the former dominates in the full R range listed). Ob- 
viously, Eq. ([l^ ) provides, strictly speaking, the refiable 
values for the kinetic exchange integral only in the limit 
U — Km ^ |tm + Kn|- All the integrals are of antifer- 
romagnetic character, since they express virtual hopping 
processes between m-neighbors, which occur (and dimin- 
ish the system energy in the second order) only when 
spins on the two sites are oriented antiparallel. The re- 
sults for the integrals ^''^ only weakly dependent 
on TV, as exemplified in Fig. |^ab, where the first two ex- 
change integrals have been shown for N = 6, 8, 10. The 



kinetic exchange is rather strong as for e.g. R = 2.65A we 
have that 42 = 0.124 eV (« 1440 K), = 1-6 meV 
(« 19 K), 41^^ = 0.44 meV (« 5 K). Such a strong 
superexchange is observe|i|Only in the two-dimensional 
superconducting cupratealj. 

The strong system-size dependence of the spin-spin 
correlation functions is shown in Fig. p^-c. This is not 
strange, since the spins correlate over much longer dis- 
tance than the range of the interactions as express collec- 
tive properties of the system (the same can be said about 
the Fermi discontinuity An^^). 

TABLE IV. The R dependence of the Fermi discontinuity 
Anfej^^^/2ij and the value obtained by the parabolic interpo- 
lation Anr". 



R/ao 


Arifcy 


A par 


2.0 


0.8264 


0.7326 


2.5 


0.6148 


0.4244 


3.0 


0.4033 


0.1582 


3.5 


0.2645 


0.0358 


4.0 


0.1791 


-0.0040 


4.5 


0.1232 


-0.0140 


5.0 


0.0857 


-0.0130 



TABLE V. The spin-spin correlation functions (Si • Si+p) 
for p = 1, 2, 3 and A'^ — 10, as a function of lattice constant. 



R/aa 


(Si ■ Si+i) 


(Si • Si-1-2) 


(Si ■ Si+3) 


2.0 


-0.2959 


0.0431 


-0.0951 


2.5 


-0.3748 


0.0845 


-0.1399 


3.0 


-0.4467 


0.1230 


-0.1793 


3.5 


-0.4968 


0.1501 


-0.2052 


4.0 


-0.5277 


0.1668 


-0.2200 


4.5 


-0.5442 


0.1754 


-0.2241 


5.0 


-0.5477 


0.1762 


-0.2175 



TABLE VI. Distance dependence of kinetic exchange inte- 
grals Jj^™' for m = 1, 2, 3, as well as the total kinetic exchange 
and the summary direct exchange integral (^ J,n). 



R/ao 


[Ry] 


r(2) 
•■'kex 

[mRy] 


r(3) 
•^kex 

[mRy] 


T{m) 
•■'kex 

[Ry] 


[mRy] 


2.0 


1.190 


37.34 


0.493 


1.227 


11.98 


2.5 


0.438 


11.67 


0.129 


0.450 


8.37 


3.0 


0.183 


3.99 


0.036 


0.187 


5.99 


3.5 


0.082 


1.39 


0.010 


0.084 


4.25 


4.0 


0.039 


0.48 


2.7- 10"3 


0.039 


2.96 


4.5 


0.019 


0.42 


2.5 ■ 10"^ 


0.019 


1.99 


5.0 


9.3- 10-3 


0.12 


4.4 ■ 10"* 


9.4 ■ 10-3 


1.27 
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V. ELECTRON-LATTICE INTERACTION, 
ZERO-POINT MOTION AND DIMERIZATION 
FROM THE FIRST PRINCIPLES 

A. General form of local electron-phonon coupling 

We can extend our |4aethod to include the local 
electron-lattice couplingpl. In the second-quantized 
Hamiltonian of general form (|l|) the positions of nuclei 
{■f^i}i=i N regarded as fixed, i.e. the ions are re- 
garded as classical objects. If their positions are subject 
to a local shifts {S^i}^^^ j^, then the Hamiltonian will 
change by the amount SH, so that now H H + SH. 
The amount of the change can be calculated as 



N 



SH 



N 



(14) 



which holds true for |(5R| <C R- Primarily, the 
nuclei shifts {(5Ri} is felt by the potential energy 
V (r) = J2i V {r — Ri) and by the Wannier functions 
{wi (r) = Wilr — Ri)}. These changes, in turn, induce 
the alteration of the microscopic parameters Eq, U, Kij, 
etc. In result, we can write the Hamiltonian change in 
the form 



SRr 



■ SRi ni + 



5U 



(<5R, - (5Rj>ta,a + ^ ^ 



^ SR^ 



{SRi — SRj) niJij . . 



SRim^mi 



(15) 



One should note that we have not included explicitly the 
more distant displacements, i.e. the terms of the type 



E 



S€a 



SRjrii, 



St 



SR, 



SRkala, 



etc. 



The last terms represent a coupling to more distant ions 
on the atomic level position, or on the nearest-neighbor 
hopping, etc. They are not included as we would like 
to determine first the derivatives S/SRi{...) from our 
above first-principle results (the subscript "0" means 
they are calculated for the periodic arrangement of ions). 
In Fig. ^ (bottom panel) we display the corresponding 
derivatives, which play the role of the local electron- 
lattice coupling constants, as a function of the inter- 
atomic distance. The scattering of the points is caused by 
the numerical differentiation. For the sake of complete- 
ness, we have plotted in the top panel the microscopic 
parameters vs. i?, for = 6 10 atoms. 



There are few unique features of these results, which 
we would like to elaborate on. First, the interaction 
terms ^ SU and ~ SK diminish the system energy when 
the system distorts. Also, dU/dR ^ —dea/dR, but ef- 
fectively \dK/dR\ overcomes dt/dR, so that the net ef- 
fects in the insulating state {hi = 1} favours the sys- 
tem distorsion (note a rather weak dependence on N). 
Furthermore, the coupling constants are relatively large. 
For example, Xa/ \SR\ = dea/dR w 0.5 Ry/ao we ob- 
tain for the distorsion \SR\/R « 0.1 (estimating the 
of zero-point motion amplitude; see below) the value 
of Aa « 2 eV for R/qq = 3 (i.e. R = 1.6l). Simi- 
larly \u = (dU/dR) \SR\ « -2.4 eV, \k ~ -1.6 eV, 
Xt ~ 0.7 eV. They represent a sizeble fraction of the 
value U — Ki « 13.8 eV. What is more important, 
Xjj/W w 0.25. However, the kinetic exchange integral 

is ~ 2.5 eV, a value close to \Xu\- This means 

that for the linear chain of hydrogen atoms the electron- 
lattice and magnetic interactions are of comparable mag- 
nitude. In the system with heavier ions the electron- 

1/2 

lattice coupling should be divided roughly by M-' , 
where Mi is the ion mass. Also, in the insulating state 
the U and K parameters (and their derivatives) should 
be roughly diminished by the relative dielectric constant 



1/2 



10, 



e of the medium. Roughly, [Mi/Mn) 
and also e ~ 5 10, so that the proportions between 
the above parameters should remain of the same order 
even though their absolute values diminish by the fac- 
tor 5 10. But this means that the analysis of the 
strongly correlated low-dimensional systems should in- 
clude on equal footing both local electron-electron and 
electron-lattice couplings. There are quite few relevant 
parameters (ta,W,U,Km,Xa,XK,Xu) in that situation, 
so the usual approach of solving Hamiltonian H + SH. by 
regarding all those quantities as free parameters of the 
model, does not look promising. Our first-principle ap- 
proach determines the value of the parameters accurately 
so their values are known for fixed R. However, needless 
to say that our method must be extended to a more re- 
alistic situation involving e.g. (Cu02)„ planar clusters to 
be applicatable to the high-Tc systems. 

The solution of the Hamiltonian incorporating the 
electron-lattice interactions requires a separate analysis 
and will not be discussed here. In the remaining part, we 
concentrate only on the evolution of the zero-point mo- 
tion and of the dimerization as a function of interatomic 
distance. 



B. Zero-point motion 

In the harmonic approximation the acoustic phonons 
for the linear chain have the dispersion relation of the 
form 



Wi. = 2 — sm — 



(16) 
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where M is the ion mass (we take the proton mass here) 
and C is the elastic constant, which can be calculated for 
longitudinal modes from the differentiation of the total 
ground state energy Eq (with inclusion of the intcrionic 
interactions), namely 



C = 



1 d^Eg 



(17) 



One should note that due to the global instability 
{dEg I dR < 0) we should place the system in a box stabi- 
lizing it (e.g. the system represents a linear ring on a sub- 
strate stabilizing its geometry). Also, the A; = is a Gold- 
stone mode, so we select the values ofA: = l,2,...,iV — 1, 
i.e. chose the center-of-mass reference system. In result, 
the contribution of zero-point motion to the system en- 
ergy is 



our rigorous analysis of the ground state properties. For 
that purpose, we define two sets of atomic basis func- 
tions A and B corresponding to even and odd lattice 
sites, shown schematically in Fig. namely 

*f (r) = *2z(r), and ^-f (r) = *2,+i(r). 

The index for each sublattice takes Nd = N/2 values 
(this is the basis halving). We construct next the two 
sets of the Bloch functions starting from this atomic ba- 
sis, which are 



N 



1/2 



exp 



i=0 



27rfcj 
1^ 



(23) 



^Ef 



k 



hujk 



Mi 



1 d^EG \ 2 
N 9i?2 



i N-1 



k=l 



In the atomic units, it takes the form 



G \ M 



1 d^Ea Y 
N dR^ 



k=l 



(18) 



7rfc\ 



To estimate the amplitude AR of zero-point motion we 
note that quasiclassically we can write for the individual 
normal mode that 



(20) 



where Ai?^ is the classical amplitude of the vibrations as- 
sociated with k-th mode. Introducing the global classical 
amplitude 



N-l 



N ^ Mujk 

k=l " 



(21) 



we obtain in atomic units that 



N 



1/2 



exp 



27rfcA 



(24) 



where the quantum number fc = 0, . . . , Nd — 1 enumer- 
ates the points in the reduced zone [Nd — N/2), and the 
normalization factors are 



■No-l 
. p=0 



AA 



COS 



2TTkp 



-I -1/2 



(25) 



The overlap integrals Sp are calculated from the pre- 
Jff^ = S (^Rff^ , where a, f3 ^ A or B, and 
j = 0, . . . , TVd - 1- Obviously, R"^^ 



scription: S°_^ = S I 



AA 



and 



Rij , the same symmetry is obeyed by all quan- 



a/3 



titles dependent on R, 

Each sublattice contains only half of the total num- 
ber of lattice sites. This circumstances leads to the 
nonorthogonality of the Bloch functions {^^^'j^) = 
SAB{k) ^ 0, even though we had before {^k^k') = Skk'- 
Therefore, the orthogonalized wave functions are 



(Ai?) 



2 a.u. 



1 



N \2m) 



1 d^EG\ ^ 



N dR^ 



-i N-l 



^ sin(7rfc/iV) 
(22) 



Approximating the summation by integration we have 
that (Ai?)^ is divergent and (Ai?)^ ~ IniV. The above 
formula expresses the first-order ^contribution to the lat- 
tice dynamics. This contributiontj appears on the top of 
the optimized energy Eg = Eg (a = amim R)- 



C. Lattice dimerization: Basis halving 

It is well known that one-dimensional systems are in- 
stable with tespect to the dimerization (the Peierls dis- 
torted phases). We incorporate the distorted phase into 



'^l=P{'^t + l*k^k. 



with 



and 



7fc 



f3k 



SAsik) 



l + v/l-|^AB(fc)P' 



\SABik)\ 



(26) 



(27) 



(l + ^l-\SAB{k)\^y 

\SAB{k)\' 



1 + Vl-|^AB(fc)|2 



-1/2 



(28) 
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The overlap integral SAB{k) for the Bloch functions can 
be related to the overlaps Sp^ between p-th neighbors 
located on different sublattices, namely 



SABik)=UkJ2 



AB |_ gBA 



■ COS 



2'Kkp 
1^ 



+ i- 



5, 



AB gBA 



„ I 2nkp 

- — sm 



(29) 



In effect, the expansion of the Wannier function in the 
basis of atomic functions has the form 



(30) 



Taking the inverse (space) Fourier transforms of the or- 
thogonalized Bloch functions we obtain the coefficients 



PfJ^j = f3p^ in the form 



dAA 



2'nkp 



BB 
P ' 



(31) 



(34) 



where the parameters are defined by the corresponding 
(primed) quantities in the atomic basis in the following 
manner 



q — r 1 



(35) 



97 



q-y .r5 
q'j^rS 



q-y,r5 
qj^rS 



(37) 



k 



Re7j. cos 



2'Kkp 



Inyyk sin 



2'Kkp 
1^ 



(32) 



BA 



N^'^^NkPk 



k 



Re7fe cos 



2'Kkp 
1^ 



luvyk sin 



2Kkp 
1^ 



Note that Pp — /3_p , in accordance with the definition 
of the relative distance Rf^^ defined in Fig. ^ 
With the help of the orthonormal basis 

{wf} a=A,B , 
i=0,...,ND-l 

we can define the system Hamiltonian with inclusion of 
all two-site interactions and all hopping properties in the 
following manner 

No-l 

= E E {''ania + UnialUial 
i=0 a=A,B 



jf3<ia 



k: 



a/3 J- tq/3 



E 



h.c. 



g7>r(5 



K' 



7^ 



O0& 
p—r 



Hr h-'p—r q—r ~ \'~^^ 

+ 4l,fif)v't]], (38) 



(33) = ^ {P'^.-'f Pl\U' + 2 ^ { 



q-y .r5 
q-y>rS 



+/5«- {Pff P'pl.J'tr] ^ 



7(5 
q—r 



(39) 



97 



g7,r(5 
(J7>r{5 



1 r ' "-^ q—r 



ol35 Y jnS 
^p—r j ^ q — r 



+ 2(3^^p^pl^ /3p-. + PfPp-,) V'tr] ■ (40) 

Those microscopic parameters are calculated first in 
the atomic basis of s-type functions regarding the size 
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of the orbitals the same on both sublattices. Af- 
ter the calculations of the electronic ground-state energy 
have been finished, we include, as before, the repulsion 
of (hydrogen) nuclei in the form (in atomic units): 

En-n - E (41) 

ia.jf3 
ia<jf3 

where each nucleus has been taken into account only 
once. Adding the electronic and inter- nuclear parts, 
we compute the ground state energy Eq/N as a func- 
tion of the average interatomic distance. In Fig. [l^ a 
we show the ground-state energy for = 2 (hydrogen 
molecule) and for the linear rings with even number of 
atoms (A^ = 4,6,8); the discussion for odd number of 
atoms requires a separate analysis, as the ground state 
configuration of nuclei is then a ionic-density wave with 
wave vector Q < tt/R. This energy is minimized with re- 
spect to both a — amin and the interatomic distance Ri 
(cf. Fig. |l^); in Fig. |l3|b we display the inverse size (a,„in) 
of the states (for A^ = 2 — 8); they are quite close to those 
obtained in the undistorted case (i.e. now Aa < 0.03), 
so that the values of the microscopic parameters such as 
U, Ki or ti are not much different in both the distorted 
and the undistorted states. Also, our analysis, while con- 
firming the existence of the Peierls distorted state in the 
finite-size systems, shows that the distorsion fades away 
with the increasing interatomic distance. 



D. Lattice dimerization vs. zero-point motion 

The zero-point motion of ions increases the system en- 
ergy and smears out their position. The dimerization di- 
minishes the system energy and, while shifting the ionic 
positions relative to each other, it leaves their locations 
sharply defined. Both effects may be pronounced in the 
finite-size systems, so the question arises how those two 
effects compete with each other. 

In Fig. |l| we compare the ground-state-energy changes 
due to the dimerization and to the zero-point motion for 
the atoms with hydrogen ionic mass. Whereas the for- 
mer decreases with the increasing size (bottom panel), 
the latter shows the opposite trend, as was discussed 
in the preceding Section. For comparison, the corre- 
sponding magnitudes of the average atomic shifts are 
displayed in Fig. |l5|. Note that the dimerization per- 
sists in the localized state of atoms and disappear only 
for R/qq > 6, when the Wannier functions can be ap- 
proximated with good accuracy by the atomic functions 
of Is type (q!~^ ~ oo). The zero-point vibrations of the 
light atoms should also contribute to the blurring of pic- 
ture of the charge densities observed e.g. in the scanning 
tunneling microscope of these nanoscopic objects. This 
type of spectroscopy should be applied to the observation 
of those effects on a local scale. 



VI. A BRIEF OVERVIEW: NEW FEATURES 

In this work we have produced a fairly complete de- 
scription of one-dimensional model system by combin- 
ing the exact diagonalization of many-fermion Hamilto- 
nian in the Fock space with the subsequent first-principle 
readjustment of the single-particle (Wannier) function. 
Electron and lattice properties have been obtained as a 
function of the lattice parameter and the microscopic pa- 
rameters have been determined explicitly. Our atjamach 
thus extends the current theoretical treatmcntsBallj to 
the strongly correlated systems within the parametrized 
(second-quantized) models by providing the determina- 
tion of those parameters (coupling constants) and, in 
turn, determining the fundamental properties of the cor- 
related state explicitly as a function of physical parame- 
ter, the lattice spacing R. Technically, we determine at 
each step the microscopic parameters taking the Wannier 
functions with an adjusted size (starting from an exact 
Wannier functions for hydrogenic-like s-states), diagonal- 
ize the Hamiltonian in the Fock space with all pair inter- 
actions and all hopping integrals included, and thus ob- 
tained ground state energy is readjusted again by chang- 
ing variationally the size of the orbitals, calculating the 
changed parameters and performing again the diagonal- 
ization in the Fock space, and so on, until the global 
minimum is reached for given interatomic distance (cf. 
Appendix C for details). This procedure is then repeated 
for each selected interatomic distance. 

Our method of approach reveals features, which cannot 
be looked into when considering only the parametrized 
models. Firstly, the atomic part of the energy is not a 
constant, as it changes widely with the changing lattice 
constant (cf. Table y). Secondly, the crossover to the 
Mott localized state of electrons, as well as the lattice- 
dimerization evolution can be studied systematically, as 
the lattice expands. Thirdly, the effective masses for both 
band and correlated states have been determined explic- 
itly. Fourthly, and probably most importantly, the in- 
clusion of all interactions (apart from the three- and the 
four-site terms) was possible, as all of them are calculated 
explicitly (otherwise, the model would contain many pa- 
rameters and become intractable or the results would be 
illegible). 

Such a procedure has been implemented so far (on a 
desktop server) for one orbital-per-site model system in- 
volving linear chains and rings containing up to A^ = 12 
atoms. However, due to the strong shrinking of the Wan- 
nier functions induced by the strong correlations, the 
results are rapidly converging with A^ and for A^ > 8 
and provide, in our view, a realistic estimate of the local 
properties of those strongly correlated systems. Obvi- 
ously, we must incorporate the screening by other than 
valence electrons, which are present in most of real sys- 
tems, as well as to extend the method to d-type orbitals 
before the analysis will become applicable to the cor- 
related 3d systems at hand. Nonetheless, our analysis 
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represents to the best of our knowledge, the first at- 
tempt to marrying consistently second- and first- and 
second-quantization aspects of the strong electronic cor- 
relations and as such should be tested in the clearest 
situatiorL. rOur method leads also to a renormalized wave 
equationOa for the single-particle wave function, but this 
feature of the method requires a separate discussion. 

The method is directly applicable to the correlated 
quantum dots, but here we avoided introducing the trap- 
ping potential, as we would like to avoid mixing phe- 
nomenological and microscopic concepts at this stage. 

In this paper we have considered only the situation 
with one-electron per atom corresponding to the half- 
fiUed-band situation in the metallic state. Therefore, the 
onset of the Mott localization washes away any Luttinger- 
liquid type of dynamics. Additionally, we have concen- 
trated on static (equilibrium) properties by determining 
the ground-state (spin-singlet) configuration and calcu- 
lating its characteristics such as statistical distribution 
function, spin-spin correlation function, the exchange in- 
tegrals, amplitude of dimerization, etc. A direct deter- 
mination of e.g. the Hubbard gap, quasiparticle prop- 
erties, pairing tendencies or spin-spin charge separation 
requires calculation of the dynamic quantities such as 
the spectral-density function and the effective interac- 
tion, particularly for the systems with one or two electron 
holes in our starting system. The evolution of these prop- 
erties as a function of interatomic distance is of funda- 
mental importance and should be carried out next. Also, 
a detailed comparison of the behavior of systems with 
odd and even number of electrons should be made. 

A separate question concerns the implementation of 
the wave-function optimization within the dynamical 
mean-field approachcj in order to obtain an explicit 
mean-field solution of a 3d model system on the exam- 
ple of the Sd Hubbard model, as a function of the inter- 
atomic distance. This would also make possible a direct 
incorporation of the band-theoretical methods with those 
emphasizing the role of local electronic correlations. We 
should see a decisive progress in this matter soon. 

Very recently, we have -calculated the effective mass 
in the interacting systemcd which corresponds to the 
quasiparticle mass in the infinite system. It is divergent 
when the calculated distribution function Uka is interpo- 
lated into a continuous distribution with a Fermi ridge. 
This results complement beautifully the results of the 
approach presented in the present paper. 
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APPENDIX A: MICROSCOPIC PARAMETERS 
IN THE ATOMIC BASIS 

In Ref.i we have expressed the parameters Ca, tij, 
U, and Kij in the atomic basis using the expansion 0. 
Here we supplement them with the formulae for the 
direct exchange integral Ji_j and so-called correlated 
hopping Vi-j. Namely, substracting ^ into the expres- 
sions Jp = Ji-j = (wiWjlVulwjWi) and Vp = Vi-j = 
{wiWi \ Vi2 \wiWj) we obtain recursively 

q qr 
q>r 

~l" PqPp—qPp—r) ^q — r ^' ^ ^ ^ PqPp—qPp—rPr { jq—r 

qr 
q>r 

+ K-r)+^T.^Pp~rJq-r' (Al) 
qr 
q>r 

= E ^^,^P-lU' + 2 E {P^qPp-q + iPpp-qPr) V'q^r 

q qr 

q>r 

+ 2 ^ PprPp-r (2 + K'^^r) , (A2) 

qr 
q>r 

where p, q and r = 0, 1,...,A^ — 1, and the expansion co- 
efficients Pp of the Wannier functions obey the symmetry 
(3p = (3_p (thus for N atoms there are N/2 or {N - l)/2 
independent coefficients when N is even or odd, respec- 
tively). The primed quantities are defined in the atomic 
basis {^j(i")}j^i ^- For example 

4- = (*.*j|l/l2 = j d''Yd''v'^*{Y)^*{Y') 

xVi2{v-Y')^,{T)^,{y'). (A3) 

These expressions are used when evaluating the micro- 
scopic parameters contained in the Hamiltonian. They 
are determined explicitly in the optimal state: Eq = 

-^inin(-R; ^ Q^min)- 
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APPENDIX B: SIGLE- PARTICLE PARAMETERS 
IN THE ATOMIC BASIS 

The expansion (||) leads to two-site terms in the atomic 
energy Ca and to and three-site terms in the hopping in- 
tegrals. For the sake of completness, we start from their 
full expressions in the Wannier basis: 



qr 
q>r 



(Bl) 



where a = Rij. The regimes for A and ^ are: 1 < A < oo, 
— 1 < /i < 1. This transformation leads to the following 
expression for 



(B8) 



where h is the ^-coordinate of the middle point of 
the Coulomb potential well caused by the k-ih ion (cf. 
Fig BnGh. Integrating with respect to ^ we obtain 



q—r' 



(B2) 



qr 
q>r 



In these expressions the atomic energy e'^ (in atomic 
units) is 



el = 



a 



2" + E 



2 / 2 

— + exp(-2aR,j) [2a + — 



(B3) 



where Rij — |Ri — Rj| and rj = |r — Rj|. Note the 
appearance of the long-range part {—2 /Rij), as one 
would have in the classical limit. 

The evaluation of t'ij is not so straightforward as one 
can see from the expression 



=To-2j2nkj, (B4) 



where = |r — Rfc|, tq represent the simple part and 
Tijk is the three-site part. The part tq is easy to calcu- 
late, since 



To 



(B5) 

The three-site part is more cumbersome, as it reduces to 
the following integral expression 



-a{ri+rj) 



(B6) 



To calculate the integral we introduce the spheroidal co- 
ordinates (A, /i) 



aX — ri + Tj , a/i 



(B7) 



1 



dAe-"°^ {[A^ (1 + {2h/af) + b/2] x 



X lo( 



{2hX/a - 1) - ^(2/iA/a- 1)2 +6 
{2hX/a + 1) - ^(2/iA/a+ 1)2 +6 



+ {3hX/a - l/2)v/(2/iA/a+ 1)2 +6 
-(3/iA/a + l/2)v/(2/iA/a - 1)2 + fej , (B9) 

where b = (A2 — 1) [l — (2/i/a)2] . This integral simplifies 
substantially if either k — i or k — j, i.e. for h = ±a/2. 
Then, taking simple limiting expression we obtain that 



dXXe' 



-aaX 



ae 



-aRi 



{aR^j+1). (BIO) 



Substituting ( BIC ) to (B4) we obtain the expression for 
the hopping integral for H2 molecule and for the_linear 
chain in the tight-binding approximation (cf. Ref.Q). 

In the general case, we have to evaluate the integral 
( |B9D numerically. For this purpose, one makes the change 
of variable A = \/t io integrate over the finite interval 
< t < 1 (one has to use a variable summation step, 
since the integrand is logarithmically divergent aXt — 1). 
When using the Simpson mcthodlla one has to evaluate 
the integrand in 300 ^ 400 points to achieve the accu- 
rancy 10~^ Ry; this procedure requires a neglegible time 
compared to that required to determine the ground state 
energy. 

Numerical calculations of all hoppings {tij}^j^^ jy 
can be accelerated by exploiting the problem symmetry. 
Among N{N — l)/2 possible hoppings for the linear-chain 
case with periodic boundary conditions {Uj = ii_j), we 
have only {N — N mod 2/A^) different integrals. In case 
of dimerization this number is doubled, as the hopping 
elements change Uj tf/^ = ^fj^^tfj^ = ^fi^- Addi- 
tionally, we use the symmetry h —h when evaluating 

'^ijk • 
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APPENDIX C: GROUND-STATE ENERGY 
EVALUATION IN THE FOCK SPACE 

For N lattice sites and the grand canonical system 
with variable number of electrons TVe = 0, 1, . . . , 2N and 
the total number of electrons with the spins up = 
0,1,..., iVe, we have the Fock (occupation- number) space 
Hi = 7ii{Ar, iVe = 0, . . . , 2iV, iV| = 0, . . . , Ne} of dimen- 
sion 4^. Additionally, if the number of electrons is fixed 
aX Nf.— N (i.e. for the case of one electron per atom) then 
the Fock space = n2{N,Ne = iV,iV| = 0, . . . ,iVe} 

has reduced dimension to dim 7^2 — ( yy^) ■ Moreover, if 
the total 2;-component of spin is S^^^ — (for N even, as 
happens in the situation then the space 7^2 {■ ■ ■}) reduces 
to the space n3{N,Ne = N,N^ = N/2}, with the cor- 
responding dimension dim Ti.^ — ^ jYy2^ • Those dimen- 
sions for = 10 are, respectively: dimTii = 1 048 567, 
dimH2 = 184 756, and dimT^s = 63 504. The implemen- 
tation of the translational symmetry does not influence 
essentially either the computing time or the memory ca- 
pacity required. 

The basis vectors in the Fock space are represented by 
the site (Wannier-state) occupancies {nia}i=i ... n as 

\v) = ...,nArt,ni|, ...,nNi) ■ (CI) 

The creation and annihilation operators are defined in a 
standard manner: 

al |..., n,,, ...) = (-I)''- (1 - n,,) \...,n,, + 1, ...) , (C2) 

Qia- \.-,nia, ...) = (-1)'^^'^ Hia- \...,nia- - 1, ...) , (C3) 

where Vicr is the number of electrons in the states preced- 
ing the i-th site (including the opposite spins if cr =i). 
We have also used the notation |. . .) =0. The above 
definitions allow for an unambigous determination of 
the matrix representation of the Hamiltonian. Simply, 
we have to determine the matrix elements {u\Ti.\v) for 
u,v € Tis- In practice, the basis vectors {\v)} are ordered 
as the combination series. Therefore, acting on each of 
them with consecutive terms of 7i, we immediatly can 
identify the result of the action as proper basis vector 
{\u)} (or zero vector). Thus, the number of operations is 
equal to the number of nonzero elements {{u\ H. 

The method allows for the construction of the Hamil- 
tonian matrix representation for which the lowest eigen- 
value is the ground state energy. The physics of the 
problem is determined usually by ~ 10% of the eigenvec- 
tors with the lowest eigenenergies. In this situationJJie 
Lanczos methodtj in the version proposed by NishinoEj is 
most appropriate. In its latter version one uses the cum- 
mulant expansion containing up to 10^ 10'^ terms for 
iV = 6 10 atoms with the lattice constant i? = 2 3ao 
(for larger R and/or smaller N the number of terms to 



be included is much smaller). The numerical accurancy 
of the result for Eq was typically 10^^ Ry. 

The final step involved the minimization of Eq with re- 
spect to a and (whenever possible) with respect to R (in 
the case of dimerization it is the shorter lattice constant 
i?i). The minimization was carried out using-j:espec- 
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make use of the translational symmetry. 
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FIG. 3. a) The space profile of the band shape for N = W 
versus R, taking into account the calculated hopping integrals 
{tm}m=i,...,5- The horizontal plane intersecting the band 
marks the Fermi level position; b) the flattening of the band 
shape with the increasing R (the dashed lines describe the 
band shape if only ti is included). 



FIG. 1. Schematic representation of the finite linear chain 
with periodic boundary conditions used in the calculations. 
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FIG. 2. Space profile of the Wannier-function evolution 
with the increasing interatomic distance. 
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FIG. 4. a) The band effective mass at the band center 
{k = 0) and b) for the Fermi wave vector k = kp, both vs. R. 
Note a general insensitivity of the results for different N. 
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FIG. 5. The profile of tlie (bare) band effective mass as a 
function of the wave vector k and the interatomic distance R 
(for N = 10). 



1.0 



0.5 



0.0 



N = 10 

a = amin 



















R/a„ = 4_l 
FVao = 3_] 








(b) 



-0.5 





Ml 



0.5 



FIG. 7. The statistical distribution function riktr for 
A'^ = 10 in the state with optimized orbitals (a) and with 
non-optimized orbitals (b); both for different interatomic dis- 
tances R. 
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FIG. 6. The Mott-criterion value ncan vs. R and for dif- 
ferent number of sites N = 6-i- 10. The horizontal line marks 
the value for bulk systems. 
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FIG. 8. Correlation functions versus distance R, charac- 
terizing the crossover from itinerant to localized state. The 
shaded areas arc drawn to illuminate the results convergence 
with the increasing R. 
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FIG. 9. dependence of the kinetic exchange integrals J^^^ 



(a) and 4;;; (b) for different N = 6 ^ 10. The results for the 
dominant integral (a) are weatly N dependent. 
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FIG. 10. Distance dependence of the spin-spin correlation 
functions (S, • Si+p) for p = 1 3 and N = 6 -i- 10. Note the 
oscillatory character with p expressing the antiferromagnetic 
correlations. 
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FIG. 11. The interaction parameters (a) and their deriva- 
tives (b) with respect to the lattice-parameter change; the 
latter express the coupling-constant strength for different lo- 
cal electron-lattice couplings. 
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FIG. 13. Ground state energy (a) and the optimal inverse 

orbital size (b) of the dimerized chain versus R, for N = 2^8 
and for the optimized both the distorsion {1 — Ri/R) and the 
orbital size q"'^. The dotted base lines mark the correspond- 
ing quantities for the H2-molecule. 




17 



0,03 



-0,02 



Zero-point motions, 




—r- 
3 



INTERATOMIC DISTANCE R/a 



0,00 



Zero - point motion : 




INTERATOIVIIC DISTANCE R/a„ 



FIG. 14. Ground-statc-cnergy changes due to the dimeriza- 
tion (bottom panel) and to the zero-point motion (top panel), 
plotted for different A'^, as a function of R. 



FIG. 15. Atomic shift due to the dimerization and the 
root-mean-square amplitude of the zero point motion vs. R 
and for different N (top and bottom parts, respectively). 




FIG. 16. Configuration of coordinates used to calculated 
the three- site terms Tikj in the hopping integral t'^j for the 
electron transfer j ^ i induced by the Coulomb potential of 
A;-th ion. 
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